Spin wave spectrum of a disordered double exchange model 
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A double exchange model with quenched disorder for conduction electrons is studied by field 
theoretical methods. By using a path integral formalism and replica techniques based on it, an 
ensemble-averaged spin wave dispersion of the localized spins is derived. It is shown that the 
spectrum of the spin wave has gaps at the multiple of the Fermi wavenumber of the conduction 
electrons in the presence of disorder, and hence, quenched disorder for electrons adds a striking 
effect to the dynamics of the localized spins. In the strong disorder limit, the present results suggest 
spin-glass like behavior due to the frustration of the exchange coupling. 



The double exchange (DE) model [J HUH has been 
attracting much renewed interest due to the discovery 
of colossal magnetoresistence of perovskite manganitcs 
|H 0- These materials show a variety of phases, in 
which spin, charge, orbital, and lattice degrees of free- 
dom and their interplay provide rich and complex behav- 
ior. Although a lot of experimental and theoretical stud- 
ies have been clarifying some aspects of complex physics 
of manganites, they still remain fascinating and incom- 
pletely understood phenomena 0- 

It is believed that the DE mechanism is essential 
to the metallic ferromagnetism of manganites. How- 
ever, several authors have pointed out the importance 
of disorder, which is inevitably included in the materials 
[3 0, El E3l ■ Especially, recent experiments have re- 
ported anomalous behavior of spin wave spectrum, which 
cannot be explained by a simple DE theory |9L 1 1 OL 1 1 lL 1 1 2| . 
Motome and Furukawa 17] have numerically studied DE 
model with quenched disorder and found that the model 
successfully describe the anomalies of spin wave disper- 
sion. They have suggested that the origin of the anoma- 
lies is the Friedel oscillation. 

Motivated by their work, we study, in this paper, the 
effects of quenched disorder on the spin wave by field 
theoretical method. Using the replica approach for the 
ensemble average over disorder and integrating out the 
electrons, we derive an effective action for the localized 
spins. The equation of motion derived from this action 
leads us to the usual one for the Heisenberg ferromagnets 
but with 2Uf oscillating terms. To clarify our ideas in 
this paper, we restrict our discussions to one spatial di- 
mension, since disorder effects are more relevant and add 
a more striking effect in lower dimensions, as claimed by 
Motome and Furukawa. 

The Hamiltonian we study is given by |l7j 

H = -tJ2 { c ]<y c i+i<y + h - c -) - Jh s i ' 8 i 
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as a classical object Sj = Srij with a unit vector rij de- 
fined by — (sin 9 j cos 4> 3 -, sin 9j sin (f>j , cos 9j ) , and ej is 
on-site disorder potential with the gaussian distribution 

P[e] = exp [- Y,j $/2g ■ 

We start with the canonical DE model without dis- 
order to have a path integral formalism convenient for 
taking disorder into account. The Hund's-rule coupling 
term can be written as Jjj $3,- s j ' Sj = {Jh S/tyCjSjCj , 
where = (c 



and 



cos 6>.j 



sin^j-e 



(2) 



The partition function is then represented by a path in- 
tegral form Z = J DcDc' e~ s , where the action is given 

by S = J Q P dr fe . c]d T c 3 -H-^Ej eje,) ' The ma " 
trix Sj in Eq. (J2J is diagonalized by a local unitary ma- 
trix Uj(r) ee S J (0 J /2,(j) J ) as U](T)SjUj(T) = a 3 . Let 
us define a fermion field Cj = UiCi in the locally rotated 
frame. For a large Jh, only up-spin fermions survive, and 
we have an effective action described solely by fermions 
Cj| ee Cj in this limit, 



S 



cj (d T + (nj\d T rij)) Cj 



-t (c^(nj\nj+i)cj+i +h.c.) - /xcjcj 



(3) 



It should be noted that in this projection the nonlocal 
transformation yields nontrivial additional terms such as 
(Ujd T Uj)^ and (UjUj +1 )^ in the canonical and hop- 
ping terms, respectively. To describe these, we have in- 
troduced a bracket notation 



H) 



TT 



cos — cos — 
2 2 

(ni\nj), 



.sinZlg-i^*-^) 



(4) 
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since this is actually equivalent to the inner product 
denotes a spin operator of the of the spin coherent state, \rij) — cos(0j/2)| |) + 
conduction electron, Sj a localized spin which we treat sin(0j/2)e l ""| j). Therefore, (rij\d T nj) = Aj(r) can be 



where St 
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interpreted as the Berry phase ^| , playing an important 
role in the dynamics of the localized spins. 

Assume that the localized spins are almost aligned to 
one direction, (n^n.,} ~ 1. In this case we introduce a 
small field variable (ni\rij) — (ni\rij) — 1, and then the 
action can be divided into two pieces S = Sq + S n , where 



Sq — 

Sn = 



J dr^2 ^jdrCj - t (cjCj+i + h.c.j - ixC-Cj 
J dT ^ \ c ] A i c i ~ 1 { c ) ( n i l n J+!) C J+ 1 + h - c -} 



(5) 

These equations are basis of our perturbative calcula- 
tions. At the leading order, an effective action W[n] for 



the localized spins, defined by e w \ n \ 
by 



(e b "), is given 



rP 

W[n] = drY] 
Jo 



9]Aj — ^Jjnj 



n 
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except for some irrelevant constant terms, where 



9j ee (c\(r) Cj (r)), J, ee t{c](r)c j+1 (r)), (7) 

and we have used the fact that | (n j | rij ) | 2 = | (1 + rij • nj) , 
namely, (rii\rij) + c.c. ~ \{rii ■ rij — 1). In the pure 
model under consideration, gj = 1 — x = go and Jj — 
It J2k>o f( £ k) cos & = Jo arc uniform constants, where 
x is the hole concentration, Sk = —2t cos k — \i is the 
transfer energy of electrons, and /(e) = l/(e /3e + 1) is 
the Fermi distribution function. It should be noted that 
in Eq. JJjJ, rij ■ rij + i ~ 1 is implicitly assumed. 

The classical equation of motion for the localized spins 
rij can be derived from the effective action Eq. JJjJl: 
The variation of the action W[n] with respect to rij, 
SW/Srij(T) — 0, leads us to 



drij 

gj —=n x {Jj_ inj -x 



Jin 



3+1/ 



(8) 



after the replacement r — > —it. In deriving this, the 
Berry phase is involved in the time derivative of the spin 
variable rij. In the pure case, this equation is just the 
one for the Heisenberg ferromagnets, giving dispersion 
relation ujk = (2Jo/go)(l — cosfe). 

Now let us take account of disorder effects on the spin 
wave spectrum. Since the disorder potential in Eq. is 
invariant under the unitary transformation Uj (r) , it gives 
rise to Hd = J2i € i c l c i m the projected action J3J). We 
now introduce m replicas to take quenched average over 
disorder. Ensemble average yields the following term 



S rf = 



® ® i a, 6—1 

m 

= ^ E E E C i+ q ,a( iuJ n)ck,b(^n>) 
a,b— 1 u n .Lj' n k,k' ,q 

xc l>-„,b( iuJ n') c k'A iuJ n), (9) 



where a,b = 1, 2, • • • , m denote replica species. Although 
the fermions in the present theory are spinless (fully- 
polarized), this inter-replica coupling in Sd as well as 
the spin-fermion coupling in S n are reminiscent of the 
theory of the spin density wave (SDW): They are ex- 
pected to play similar roles to the Coulomb interaction 
and phonons in SDW, respectively. The wisdom obtained 
so far tells that 2kp oscillating mode of phonons makes 
a gap in the fermion sector, where Uf denotes the Fermi 
wavenumber. On the analogy of this, we expect that the 
2Uf mode of localized spins in S n would be potentially in- 
volved with a gap formation of fermions, although we will 
treat S n separately from So + Sd, in what follows. Based 
on these observations, let us first examine the mean field 
ground state of the fermions within the action So + Sd, 
and next compute the coupling constants of the localized 
spins, gj and Jj in Eq. Q. We have to mention here 
that Sd is nonlocal with respect to r, and therefore care- 
ful treatment of the inter-replica interaction is needed. 
Namely, we divide the summation over the Matsubara 
frequencies into two parts: 



+ f E E E E C fe+ g ,a(* W «i +iVn)ck,a(iu ni )cl,_ q j(iui n2 - iUJ n )c k ^ b {iuj n2 ) . (10) 
a, b Wn^Ow^jWnj q,k,k' 

Here, uj n — 27rn//3 is bosonic Matsubara frequency. It should be noted that a = b terms of the first line of the right 
hand side above vanish identically due to the Fermi statistics. Introducing two kinds of auxiliary fields leads us to 
the action which is bilinear in the Fermi fields, 

Sd = r E E E Q?,ta( iw n)Q-?,at(- !w ») + 7T E ^ Pq( iuJ n)P-q(~iUJn) 
" u>„ q a-^b " w„/0 q 
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+^ ^ ^ ^ \Q q .}, a {ibj n )c\ +qa {iuj ni + iu„)c kib (iu n J + Q- q , ab (-iu n )cl_ qb (iuj ni - iu)n)ck,a(iu ni ) 

Ul n «? k a^b 

+ ~2 [ P «( iWn )4+g, Q ( iw ni +^ .(11) 

u) n =jtG q u> ni k a 

I 



Let us now assume a static (iu> n = 0) saddle point solu- 
tion with the following replica symmetric form 

Q q ,ab(iu n ) = 2iA(35 n , 5 q ,2k F , (12) 

whereas P q {iui n ) = since iui n = is prohibited in the 
summation in Eq. This form of the mean field so- 

lution is similar to the one in the replica theory for spin 
glass models 0] . Although we have mentioned the simi- 
larity to the theory of SDW, the present replica coupling 
is repulsive. Therefore, we cannot expect a stable mean 
field ground state if the number of species m is an inte- 
gers m > 1: This is the reason we assume an imaginary 
Q in Eq. lO- Substituting Eq. (O into Eq. |TTJ) , we 
have 

s d = - A [ c i+2k F A iuj ™) c ^ iuj n) + h - c 

a^b tJn k 

2m(m — 1) . o . 

i -A . 13 
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Due to the assumption of imaginary Q, the A 2 potential 
above has actually a wrong sign. Nevertheless, in the 
replica limit m — > 0, we have a mean field solution, as we 
shall see momentarily. 

Now let us determine the parameter A by minimizing 
the effective action for A. To this end, we decompose the 
Fourier modes into several pieces, 



Cj,a(iu n ) = e lkj c ky a(iuj n ) + ^2 u' kij c kia (itA} n ),(li) 

\k\<k F -A \k\<A 

\k]>k F +A 



Where 4,o = ( C l+k F ,a' C l-k F ,a) and U k,j = 

( e i(k+k F )j ^ e i(k~k F )jy Then the fermion action becomes 



s + s d = 



a.b 



y^ c\ M (iu> n ) (iU) n - Ek) Cfc,&(lW„) - c\ a (iLO n ) [D k (iui n )Sab + M(l - S a b)] Ck,b(i0Jn) , 

(15) 



\k\<k F -A 
fc|>fc F +A 



fc|<A 



where 



D k (iuj n ) 

M : 



iuj n — VFk 

iuj n + vpk 



A 
A 



(16) 



Here we have linearized the transfer energy of electrons 
as Sk±k F ~ ivpk with vf = 2tsinfc^. Note that the 
matrix Dk5 ab + M{\ — S ab ) can be diagonalized in the 
replica space by a global unitary transformation, giving 
(m — 1) eigenvalues Dk — M and 1 eigenvalue Dk — (1 — 
m)M. Therefore, integrating out the Fermi fields yields 
/X>c2?c t e-( So+ ' Sd ), where 

= - 2m(m ~ 1)0 A 2 - (m - 1) InDetQP - M) 



-F(A) 



F(A) 



InDet (D- (1- m)M) . 



(17) 



It is readily seen that Dk(iu> n ) — M has eigenvalues 



iw„ ± Bfe(A) with £*(A) = ^(v F k) 2 + A 2 . Then, the 
variation with respect to A gives 



dF J 4(1 m) 

<9A = | 5 



■In 



Vf 



+ 0(m) W18) 



where we have assumed that the temperature of the sys- 
tem is much lower than A, i.e, /3A <C 1. This equation 
gives no solution in the normal case with m > 1 , since the 
inter-replica interaction is repulsive, as mentioned above. 
However, it turns out that we obtain non trivial stable 
solution in the replica limit m — > 0, 



A = 2v F A- 



2v F Ae~ 



(19) 
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This is due to the change of the sign of the A 2 potential 
at m = 1. 

Now we calculate the effective action Eq. (0 in 
the presence of disorder by using ensemble-averaged 
correlation functions in Eq. J7J. To this end, we 

need (cj(T)cj(Y)} = lim m _>o(c] i( T ) c j,i( T )) m the replica 
approach. In order to calculate this in the replica- 
diagonal frame introduced previously, we have to note 
that the above correlation function is equivalent to 

lim m ^ ^ S Q ( c la( T ) c J.«( r ))' and therefore, invariant 
under the global transformation. As we mentioned al- 
ready, this transformation yields m— 1 eigenvalues D+M 
and 1 eigenvalue D + (1 — m)M, the latter of which ap- 
proaches the former in the replica limit m — > 0, and 
therefore, this observation enables us to keep just one 
replica species, say, a = 1 with D + M to calculate the 
above correlation function. To simplify the notation, we 
will suppress the replica indices hereafter. 
Based on 



CO 



E 



JHj-j'). 



\k\<k F —A 
fc|>Jc F +A 



tuj n - e k 



^ tTu k jml tj (c k (ioj n )cl(ioj n )) 

|fc|<A 



. (20) 



and assume that A is much smaller than the temperature 
under consideration, we end up with 

9i = 9o - g osc cos(2k F j), 

Jj = J Q -5J-J osc cos[k F (2j + l)}, (21) 



where each term can be computed as 

A _ A 
E k nv F 

5 J = 2tsmk F sin A; | 1 — — I , 



„ A A . , f Av F \ 

= 2 > — — = arcsmh — — 

^ E k irv F V A J 

0<fe<A 



0<fc<A 



josc = 2t ^ cosfc 



0<fc<A 



A 

E~ k 



(22) 



Thus we have derived a classical ferromagnetic spin- wave 
action © with Eqs. (|21|) and (|22|) . which yields the equa- 
tion of motion JSJ. Eq. (|21() tells us that nonmagnetic 
impurities of conduction electrons give rise to a 2k F os- 
cillating component to the exchange coupling of the lo- 
calized spins. This results in the gap opening in the spin 
wave dispersion io q at q — hp, 2k F , ■ ■ ■ mod(7r), which are 
controlled by the Fermi wavenumber k F of conduction 
electrons. We show in Fig. ^ a schematic spin wave 
spectrum for a x = 1/3 doping case. As expected from 
the exchange coupling with the 2k F oscillation, there ap- 
pears a gap at q — (2/3)7r. In addition, the dispersion 
relation has a small gap at q — 2k F — (l/3)ir mod tt. In 




2k, 



k F 



FIG. 1: Schematic illustration of the spin wave dispersion for 
the x = 1/3 concentration case in the periodic zone scheme. 
Gaps open at q — k F = 2/3n and 2k F = 1/37T (mod n). 



general, when the Fermi wavenumber is commensurate, 
k F = (m/n)ir, there appears n gaps in the spin wave 
dispersion. It should be noted that 2kp oscillation is in- 
cluded not only in the exchange coupling but also in the 
Berry phase term, both of which determine the size of 
gaps. As far as we have studied the equation of motion 
numerically, gj tends to make the gap at q — kp larger 
and the other gaps smaller. Details will be published 
elsewhere. 

In a strong disorder limit, we expect relatively large A 
of order A. Although the present formalism is a weak dis- 
order approach, it may be suggestive to consider the case 
A — > oo. Setting A ~ kp, we have gj ~ 2(1— x) sm 2 (kpj) 
and Jj ~ — s'm[kp(j + 1)] sin(kpj). These equations tell 
that the exchange coupling could be frustrated, which 
may be responsible for the glass phase of the manganites 
observed actually by experiments 20]. However, in this 
case, quantum treatments of the localized spins would be 
required. 

So far we have studied a spin wave dispersion in one 
dimension, which has gaps for any disorder strength g. In 
higher dimensions, there exists a critical value g c (> 0), 
below which the dispersion has no gaps. Details will be 
published elsewhere. 
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